function data = load_Rb_D2Line()
    try
        data =load([PUMPING_PROJ_PATH '/phys/+CellSystem/@VacuumCell/VacuumCellRbD2.mat'], 'res_absorption', 'freqList');
        absorption = data.res_absorption;
        freqList = data.freqList;
    catch
        [absorption, param] = CellSystem.VacuumCell.vacuum_cell_absorption_RbD2();
        freqList = param.Frequency;
    end
    plot(freqList, absorption(1, :), 'r.-', freqList, absorption(2, :),  'b.-', freqList, sum(absorption, 1), 'k.-')
end
